Temperature-dependent quantum electron transport in 2D point contact 
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We consider a transmission of electrons through a two-dimensional ballistic point contact in the 
low-conductance regime below the 0.7-anomaly. The scattering of electrons by Friedel oscillations 
of charge density results in a contribution to the conductance proportional to the temperature. The 
sign of this linear term depends on the range of the electron-electron interaction and appears to be 
negative for the relevant experimental parameters. 
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I. INTRODUCTION 

In recent years, the effects of electron-electron inter- 
action on the conductance of low-dimensional ballistic 
contacts have attracted a considerable interest. This 
is mainly due to the attempts to explain the conduc- 
tance plateaus at 0.7(2e 2 //i) in quantum point contacts^, 
and 0.5(2e 2 //i) in quantum wires- Usually these effects 
are explained by an existence of a localized state in the 
contact 3,4 . The mechanism of localization is not fully 
understood, and different scenarios of formation of such 
a stated— were proposed. One of them involves electron 
backscattering from the oscillations of the electron den- 
sity in quantum point contacts^. This mechanism is in- 
herent to any constriction, and the oscillations of electron 
density are actually observed in experiments^. 

The formation of the plateaus in the gate-voltage de- 
pendence of the conductance is not the only possible ef- 
fect of electron scattering by the Friedel oscillations in 
2D systems. Recently, it was shown that scattering by 
Friedel oscillations in a two-dimensional conductor with 
impurities results in a strong temperature dependence of 
the conductivit y 10 ' 11 . It was also predicted that this scat- 
tering may give rise to a zero-bias anomaly of tunneling 
into the edge of a 2D electron gas^. 

Most of theoretical papers^ - — dealt with interaction 
effects in the narrowest part of quantum point con- 
tacts (Fig. Q} by considering them as ID channels and 
took into account only a few lowest transverse quantum 
modes. Therefore it is not clear how the transition to 
the continuum of quantum modes in the electrodes takes 
place and whether the interaction effects in the transition 
region outside the constriction play a role. Meanwhile it 
is well known that Friedel oscillations in a 2D electron gas 
fall down with distance x from a planar barrier accord- 
ing to the law a; -3 / 2 , i.e. they penetrate deep into the 
electrodes. Hence their contribution to the conductance 
may be significant. 

Recently, we considered contacts much wider than 
the Fermi wavelength and obtained the temperature- 
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dependent contribution to the conductance**, positive 
weak-field magnetoresistance^, and the shot noise^ due 
to electron-electron interaction in the semiclassical ap- 
proximation. The theoretically predicted G(T) and 



constriction 



FIG. 1: Sketch of a realistic quantum point contact. The blue 
triangles show the gates that form the constriction, and the 
red arrows show the direction of motion of electrons. 



G(H) dependencies^ 3 .^ are in a good agreement with 
the experiment s 16 : 17 . However the linear in temperature 
contribution to the conductance persists even for con- 
tacts of the width approaching the Fermi wavelength 17 , 
where the scattering by quantum oscillations of electron 
density should be considerable. Therefore it is of inter- 
est to compare the contribution from this quantum effect 
with the semiclassical one. To this end, we calculate the 
temperature-dependent correction to the conductance for 
the case of small contact size a -C and consider the 
interplay between the single-slit diffraction and the inter- 
action. As we assume the interaction to be weak, the for- 
mation of a localized state in the contact is irrelevant to 
our problem. We do not address here the physics related 
with the 0.7-anomaly and focus on the low-conductance 
regime. 

The electron scattering by the Friedel oscillations re- 
sults in a cusp in the probability of transmission through 
the contact at the Fermi surface, which leads to a lin- 
ear temperature dependence of the conductance similarly 
to the semiclassical case. However the sign of this lin- 
ear correction depends on the the competition between 
a negative contribution from the direct interaction and a 
positive contribution from the exchange interaction. The 
dependence of the absolute and relative corrections to the 
conductance on the contact size is also different from the 
semiclassical one. 

The paper is organized as follows. In Section II, we 
present the model and describe our general formalism. 
Section III addresses the case of noninteracting electrons, 
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and Section IV describes the perturbation theory. Sec- 
tions V and VI present the results for a point-like inter- 
action and a generalization for an interaction of a finite 
range, and Section VII contains the discussion of the re- 
sults. 



II. GENERAL APPROACH 



We consider the effects of electron-electron interaction 
on the conductance of a narrow short contact at non-zero 
temperature. We assume that electron-electron interac- 
tion is weak so that it can be treated perturbatively. 

As the Friedel oscillations die out at a large distance 
vp/T from an obstacle, we are mainly interested in scat- 
tering processes that occur in the regions outside the con- 
tact in the leads and do not focus on the exact dynam- 
ics of an electron in the narrowest part of the constric- 
tion. Therefore we consider an extremely short contact, 
namely, we use typical single-slit diffraction model - a 
gap of width 2a <C Xf in an one-dimensional barrier sep- 
arating two half-planes of 2DEG (Fig. d]) 18 . This model 
geometry allows us to avoid dealing with an infinite num- 
ber of discrete transverse modes and to use instead the 
continuous representation. 

Note that Friedel oscillations far from the barrier do 
not depend on the exact shape of the confinement poten- 
tial because they are formed by electrons near the Fermi 
level with almost normal incidence on the barrier. Hence 
a smooth barrier potential should result only in a shift 
of their phase, which would not essentially change the 
correction to the conductance (see Appendix B). 

We obtain the conductance by using a classical Lan- 
dauer approach 1 - and write the conductance as a sum of 
transmission coefficients 



°-T/i(-fO£w-< 



(i) 



Here g s is a spin degeneracy and f(k, q) is the transmis- 
sion amplitude from mode with the wave- vector k in the 
left half-plane to the mode with the wave- vector q to the 
right half-plane. 

First of all we calculate t = to and G — Gq for nonin- 
teracting electrons. A weak electron-electron interaction 
results in a scattering of electrons by the Friedel oscilla- 
tions caused by the contact boundaries. We consider the 
oscillations arising from the barrier as one-dimensional 
and neglect their distortion by the gap because this ef- 
fect is of higher order in the contact size. The incident 
electron is scattered by the Friedel oscillations before 
and after passing through the contact, which results in a 
correction to the transmission coefficient of the contact 
t(k, q) = to(k, q) + 5t(k, q). The correction to the trans- 
mission coefficient may be obtained by expanding the 
perturbation of the wave function Sip in plane waves. To 
calculate this perturbation, we solve a Schrodinger-type 
equation 




FIG. 2: The model of quantum point contact used in our 
calculations. Red line shows a process of electron scattering 
by the Friedel oscillations that affects the conductance. 



ip(r) = e ip(r) 



(2) 



in an iterative way-22, and obtain Sip(r) in the lowest order 
in the interaction. 

Now we discuss the interaction potential induced by 
the Friedel oscillations. It is a sum of a direct term and 
an exchange one^i V e f / (r) = Vh (r) — Vf (r) , where 



Vff(r)=g s J dr x U ee (r -rx)n(rx,rx), 
V F (r)ip(r) = f driU ee (r -n)n(r, n)ip(ri), 



(3) 



(4) 



where n(r,ri) = (ip + (ri) ip(r)), ip + and ip are electron 
creation and annihilation operators, and U ee (r — ri) is 
the potential of the electron-electron interaction. Typi- 
cally, it is the Coulomb interaction screened by the two- 
dimensional electrons and by the gate. The coefficient of 
spin degeneracy g s appears only in the direct term be- 
cause it involves interactions between electrons with both 
spin directions while the exchange interaction is possible 
only for electrons with the same spin. 

Equations (|3]) and Q result in a correction to the wave 
function in the form Sip — Sipu — Sipp and, accordingly, 
in a correction to the conductance 



SG = SGh — 5Gf- 



(5) 



The negative sign in the exchange term is explicitly 
shown here. 



III. NON-INTERACTING ELECTRONS 

In the absence of electron-electron interaction the con- 
ductance calculation reduces to the standard problem of 
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diffraction by a narrow gap of width a <C Xf ■ For the 
three-dimensional case it was considered many times^ 
and the conductance of a small three-dimensional ballis- 
tic contact^ was found to be proportional to the sixth 
power of the contact size G oc (kpa) 6 . 

To the best of our knowledge, the two-dimensional 
problem in the limit of a <C Xp was considered only 
once for a specifically designed model of the contact 2 ^. 
In the limit of fcj?o <C 1, these authors obtained G oc 
1/ ln 2 (/cp-a), which is unphysical. Therefore we recalcu- 
late this quantity using the solution of the problem of 
diffraction from a narrow slit obtained many decades ago 
in optics^. 

We use an approach similar to Sommerfeld— and re- 
duce the solution of the Schrodingcr equation ([2]) with 
V e ff = to a boundary- value problem. The total wave 
function may be presented in the form 




^o(r')+^(r'), 



< 
>0, 



(G) 



where ipo is the wave function in the absence of the gap 
and ipt is the lowest-order correction in the gap size. The 
zero-order wave function %po obeys zero boundary condi- 
tions both at the barrier and the gap, while the correction 
ipt obeys the zero boundary condition at the barrier and 
a nonzero boundary condition at the gap 



(v 2 + k 2 )Mx,y) = o 
V't(^,y)U=o,j/e(-o,a) = x(o,y), 



(7) 



where k 2 = 2ms. 

We assume that the incoming plane wave "0i( r ) = 
^Jmjk x e *(*»a5+feww) w ith k 2 + k 2 — k 2 falls on the contact 
from the left in the x direction and find the boundary 
condition x{v) self-consistently using the continuity of 
the derivative of the total wave function at the gap (see 
Appendix A). Expanding the transmitted wave function 
in plane waves allows us to obtain transmission coefficient 
for noninteracting electrons 



-ya 2 \Jk x q x 



i (k, q 

We substitute it in Eq.([T]) and obtain 
Go 
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y H 128 F 
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El 



(8) 



(9) 



This contact-size dependence is more physically plausi- 
ble than that of Ref. [24] because it corresponds to the 
two-dimensional analog of the Rayleigh scattering of light 
by small particles 22 . Indeed, the conductance is propor- 
tional to the square of the two-dimensional particle vol- 



IV. PERTURBATION THEORY 

Now we take into account a weak electron-electron in- 
teraction. This interaction leads to a scattering of elec- 



trons by the Friedel oscillations induced by the barrier 
and results in a correction 5t(k, q) to the transmission 
coefficient. We substitute it in Landauer formula ([!]) and 
obtain the correction to the conductance in the case of 
weak interaction 



SG 



els 
2^ 



ds 



^|t (k,q)|Im«(k,q). (10) 



Here we take into account the fact that to is an imag- 
inary quantity. The correction to the transmission coeffi- 
cient is conveniently expressed in terms of the correction 
to the wave function by expanding it in plane waves. 
The wave function is found by solving the Schrodinger 
equation ([2]) in the lowest order in the interaction. To 
this end, we isolate the term with V e ff in the right-hand 
side and substitute the unperturbed wave function into 
it. The solution is given by 



<^(r)= / dr' g(r,r')V ef f(r')i>(r'). 



(11) 



Here V e ff(r') is the scattering potential produced by 
the Friedel oscillations, g(r, r') and ip(v') are the single- 
electron Green function and the total wave function for 
noninteracting electrons. We assume that the electrons 
are incident on the contact from the left and we measure 
the total current on the right, where x > 0. We are in- 
terested in the entire range of values of x' e (— oo, oo) 
because we consider the scattering by the Friedel oscil- 
lations on both sides of the contact. Similarly to the 
wave function ([S]), the one-electron Green function may 
be written in the form 



#(r,r') = g t {r,r'), 



x 1 < 



g(r,r')=5o(r,r')+^(r,r / ), x' > 0, 



(12) 



where go is the Green function in the absence of the gap, 
and gt oc a 2 is the second-order correction in the gap 
size. We calculate gt similarly to ipt (see Appendix A) by 
solving the system 



— (V 2 + fc 2 ) 3 (r,r')=<5(r-r') 
Zm 

g(r,r')\ x=0 ,ye(-a,a) = x(0, y, x', y'). 



(13) 



We substitute © and (TT2")) into (fTTj) and obtain in the 
lowest order in the contact size 

^(r)= / dr'g t (r,v')V eff (r')iPo(r') 



x'<0 

+ J dr'g Q (r,r')V eff (r')Mr')- (14) 

x'>0 

The first term corresponds to electron scattering by 
Friedel oscillations in front of the contact, and the second 
one - behind it. 
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We substitute the expressions for the interaction po- 
tential in the Hartree-Fock approximation §3§ and Q 
into (fl~4"f and obtain the conductance as a sum of direct 
and exchange terms 5G = SGh — SGf. Then we substi- 
tute ipo, ipti 3o> and gt into the resulting expression and 
after some simplifications obtain the conductance for an 
arbitrary interaction potential in the form 



oo oo 



SGh = ~9^^~T??Tr a4 ( d£ ( — -J- ) k 2 I dyi I dx\ 



h h 2 16 



de 



x J dx'U ee (x' - x 1 ,-y 1 )n(x 1 ) J dq y sin(2q x x') 

-oo 

(15) 



6G f 



h h 2 



de 



de 



dyi 



n ST oae (e) 










6 





FIG. 3: Cusp in the transmission coefficient. 

V. POINT-LIKE INTERACTION POTENTIAL 

Consider now the case of a point-like interaction poten- 
tial U ee (x' — x\, —yi) — Up 5{x' — x\) S(yi). A comparison 
of Eqs. ([13) and flT5J| shows that SG H = g s SG F . There- 
fore SG — (g s — 1) SGf. Upon an integration over xi, y\ 
and q y , one obtains the correction in the form 



x / dxi J dx'U ee (x' -xi,-yi)n(r',xi,yi + y') 
o o 

oo 

x / dq y sinf^x') cos(q y yi) cos(q x Xi) (16) 



We use the coordinate transform y\ — > y\ + y 1 to make 
the interaction potential independent of y' and then inte- 
grate over y' . This transform results in the independence 
of the Friedel oscillations of density on y 1 because we ob- 
tain the corrections in the lowest approximation in the 
contact size and use the unperturbed wave functions in 
the absence of the gap to calculate n (see Appendix B). 



n(r',x 1 ,y 1 +y') 



1 

2^ 



, , df(p)\ 



Ji(pV( x ' - xi) 2 + vl) Ji(pV(x' + xx) 2 +yf) 



y/(x' -xi) 2 +yf 



^{x' + xi) 2 +y\ 



(17) 



By setting x' — x\ and y\ = in this expression, it 
is easy to obtain the electron density n(xi) — n(ri,ri) , 
which is responsible for the direct interaction term and 
depends only on one coordinate 



h 2 1 

n{xi) = — 

47r 2ir 



P^^- (18) 
2xi 



5G=[l-g s ]g s - 1 - 2 -aU p 



de 



M 

de 



dx n{x')Ji(2kx'). (19) 



With ro(x') from (|18[) substituted into this expression, it 
is easily seen that the main contribution to it is given by 
values x' ~ v F /T, i. e. by the "tail" of the Friedel oscil- 
lations far from the barrier. We integrate over x' and p 
and calculate the total transmission coefficient ST(e). It 
is a sum of two parts 5T const (e) + ST osc (e) formed by the 
constant and the oscillating part of the electron density 
(fl~8"|). respectively. The term ST const (e) is a smooth func- 
tion without singularities, whereas ST osc (e) has a cusp at 
the Fermi surface of the form 



ST n . 



,(e) oc — 



T 

Ef~ 



1 



(20) 



At T/Ef <C 1, the derivative of the last term with respect 
to e/Ep tends to 2 at e = Ep — and to 1 at e = Ep + 0. 
This cusp (see Fig. [3]) results in a linear temperature 
dependence of conductance 



SG = 



g s \ 9s 



128 



h h 2 



T 

Bp"' 



(21) 



Alternatively, this temperature dependence may be at- 
tributed to the temperature-dependent cut-off length of 
the Friedel oscillations (see Fig. 1). 



The first term here presents a uniform charge density, 
and the second one describes its oscillations with a pe- 
riod (2kp)~ 1 at large distances from the barrier that de- 
cay according to the law x -3 / 2 at zero temperature. At 
nonzero temperature, they exponentially decay at a char- 
acteristic length vp/T. 



VI. ARBITRARY INTERACTION POTENTIAL 

A. Direct interaction 

An isotropic finite-range interaction potential is conve- 
niently described by its Fourier components U p (p), which 
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depend only on the absolute value of p. Similarly to the 
case of a point interaction potential, the substitution of 
the two terms in (fT8| that correspond to the constant 
and the oscillating parts of the charge density into ([15)) 
results in a sum of two terms SGh = SG h, const + SGh,osc- 
The first term is easily calculated and equals 



The first term here presents the contribution linear in 
temperature and is proportional to the Fourier compo- 
nent of the interaction potential dX2k F while the second 
one presents the temperature-independent contribution 
and vanishes if U p (p) is a constant. 



+ o(e- E "l T \ (22) 



Here U P (Q) is the Fourier transform of the interaction 
potential. After some simplifications, the second term 
may be brought to the form 



B. Exchange interaction 

The substitution of the two terms of (flTl) into (TTB]) 
gives the exchange contribution to the conductance in a 
form 5G F = SG FtCO nst + SG FtOSC in analogy with SGh- 
The first term is easily calculated and equals 



SG 



e m 1 



H.t 



H fr 2 128vr 



dk 



dk 



OO 

o 

where we have introduced the notation 
(-pi + 2k) 



(23) 



h(p>k) = 



Akp 



Ap 2 -p\ 



Pi 



y/Ak 2 



pi 



- Pl +2p)9( Pl -2k) k Up 2 



Pi 



pi + ^ P 2 - 4k 2 p V pi - 4fc2 ' 



(24) 



h(p<k) = 



(-pi + 2p) 
Akp 

2p) 6(- Pl - 



Ap 2 
2k) 



pi 



Pi 



y/Ik 2 



pi 



ppi 



ksjAk 2 - p\ pi + ^Jp\ - Ap 2 



(25) 



The quantity I\(k,p,pi) has singularities at p\ — 2k and 
Pi = 2p, and the derivatives of the distribution function 
in (|23[) cut out narrow intervals of k and p of width ~ 
kpT/Ep near kp. As we assume the potential U p to be a 
smooth function of p at the scale T /vf, we can isolate the 
singular part of the integrand and substitute U p {pi) — 
U p (2kp) in it, while setting Ii(k,p,pi) = I\{k — kp,p = 
kF,Pi) in its regular part 

U p (pi) h(k,p,pi) « U p (2k F ) h(k,p,pi) + 
+ [U p ( Pl ) - U p (2k F )] h{k = k F ,p = k F , Pl ). (26) 

We calculate both terms and obtain the correction due 
to the direct interaction in the form 

2k F 

,2 ^ ™ 1 ,4 4 f . U p ( Pl )-U p (0) 



e m 1 



SG F 'const — ~9s^~ To — kpa^ 



2k r 



x J dpi U p (pi)pi axccos 



Pi 
2k f 



O 



L-**/t\ 



(28) 



After a simple rearrangement, the second term may be 
brought to the form 



e m 1 



<5Gi?ose — 9s — To. 

OO 

9f 



x dp 



dp 



h h 2 32tt 2 

p fc 

dp y I dq y Up(p y -q y ) In 



q x +Px 



q x ~Px 



(29) 



This term has a singularity at p x — q x because the 
backscattering of electrons is most efficient if the x com- 
ponent of the electron momentum q x coincides with the 
wave vector p x of the Friedel oscillations. We write the 
integrand as a sum of two terms, one of which has a 
singularity at p x = q x and the second one is a regular 
function, so that one may set k = p = k F in it to obtain 



U p (p y - q y ) In 



qx +Px 



q x ~Px 
+ [U P (p y - q v ) - U p (0)} In 



U p {0) In 

q x + Px 



Qx +Px 



q x -p x 



q x —p x 



p—k—kp 



(30) 



We perform the integration in (|2^|) , sum the result with 
(f28f and obtain the correction due to an exchange inter- 
action in the form 



T 
E~f 



where 5G Fi t=o is a temperature- independent quantity 
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given by an integral 



2k r . 



em 1 



h K 2 32 



^—kpa 4 " J dpi U ee {pi) 




Pi I Pi 

7r — — arccos — — 
2k F \2k F 



(32) 



where K and E are full elliptic integrals of the first 
and second kind. The temperature-dependent correc- 
tion to the conductance in (l3lj) is determined by the 
long-wavelength component of the interaction potential, 
which is typical for the exchange interactio n 10 ' 27 . How- 
ever 8Gf,t=o is determined by all components of U p from 
to 2kp. Long- wavelength components contribute to 
(|3"2")l with a positive sign and short-wavelength compo- 
nents contribute to it with negative sign, so that the in- 
tegral is zero if U p is constant. 



VII. DISCUSSION 

The summation of © , (|2"T1) and (|3"Tj) gives the full con- 
ductance G = Go + SGh — SGf in the form 



G 



9 

e 7r 
5s ¥l28" Y 



5Gt=o 



e 2 m 1 T 
+9s J¥l28 kpai [Up{0) ~ 9sU p {^f)\ (33) 

Here SGt=o is the temperature-independent contribution 
that results from a small interaction-induced change in 
the Fermi level. To calculate the relative correction from 
the interaction, 8Gr — SGf should be divided by © (we 
are interested only in the temperature-dependent term 
and g s =2) to give 

SG T T 

—^ = v 2 [U p (0) - 2U p (2k F )} — , (34) 

Lro tip 

where v 2 — mj (irti 2 ) is the spinful density of states. The 
relative correction linearly depends on temperature and 
is much larger than the standard Fermi-liquid T 2 correc- 
tions. It is a consequence of the cusp in a transmission 
coefficient (Fig. [3]). This is essentially the same temper- 
ature dependence that was obtained previously for the 
correction to the conductance of wide ballistic contacts 
and was due to the collisions of the injected and incident 
electrons in the leads^. 

The absolute value of the temperature-dependent cor- 
rection to the conductance is proportional to the fourth 
power of the contact size 5G oc (fcj?a) 4 and the relative 
one 5G/G does not depend on a. This is in contrast 



with the correction for wide contacts SG S 



which is 




0.2 0.4 0.6 0.8 1 



FIG. 4: The sign of the correction resulting from Friedel os- 
cillations vs. kpd and K-^d. The grey region corresponds to 
positive sign, and the blue, to negative. 



cutoff length due to a scattering by impurities or a finite 
size of the sample, so that SG/Go\ sem i is roughly pro- 
portional to Go- Should this correction be extrapolated 
to narrow contacts, it would be proportional to G 2 be- 
cause both the number of injected and incident electrons 
is proportional to Go . Hence the correction from Friedel 
oscillations must dominate at small contact widths. 

The quantum correction from the Friedel oscillations 
is more sensitive to the shape of interaction potential 
than the semiclassical one. In particular, its sign is deter- 
mined by the factor [U p (0) — 2U p (2kp)] , which is due to a 
competition between a positive contribution from the ex- 
change interaction and a negative contribution from the 
direct interaction. This factor arises in several theories 
of scattering by Friedel oscillation o 8 ' 10 ' 11 and is positive 
for long-range and negative for short-range interactions. 

Consider the most typical case of the Coulomb poten- 
tial statically screened by a metallic gate parallel to the 
2D electron gas and the electrons in the gas itself. If the 
distance between the gate is d and the dielectric constant 
of environment is Sd, the interaction potential between 
the electrons in the gas is given by (see Appendix C) 



U(q) 



Aire 2 



e d {coth\qd\ + l)\q\ +47re 2 ^ 2 ' 
which leads to the correction of the form 



(35) 



SGr 
Go 



T 
~E~f 



2nid 
1 + 2n 2 d 



2k 2 



k F [coth(2k F d) + 1] + n 2 



, (36) 



where k 2 is the inverse screening length. Figure 4 shows 
the regions in the (kpd, n 2 d) plane where the correction 
is positive or negative. 

Unfortunately, we are unaware of detailed measure- 
ments of the temperature dependence of the contact con- 
ductance near pinch-off. The estimation for realistic val- 
ues of distance to the gate d — 100 nm, the inverse screen- 
ing length k 2 = 2ire 2 v 2 /Ed = 1.93 x I0 6 cm~ 1 and elec- 



proportional to (kpa) ln(/ c /a), where l c 3> a is a large tron density 



10 



-10 



cm 2 for T 



IK results in 
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5Gt/Gq = —15%. Measurements of conductance in the 
region of the 0.7-anomaly, indeed, reveal a negative slope 
of G(T) dependence 28 . However a quantitative compari- 
son with our predictions is not possible. 

In the multichannel regime, this negative correction 
should be suppressed by the positive one from the scat- 
tering of oppositely moving electrons^, and one should 
observe a change in the sign of the slope of its temper- 
ature dependence. However if the actual ft 2 is smaller 
due to the low electron concentration, the sign of the 
correction from the electron-electron interaction may re- 
main positive for all contact sizes. Therefore one may 
estimate the actual screening length from the sign and 
slope of the temperature dependence of conductance for 
narrow contacts. 

Though our results were obtained for a sharp potential 
in a form of an infinitely narrow and very high barrier, 
they should also survive for reasonably smooth poten- 
tials. It is easily seen that the smoothness of potential 
results only in a phase shift of the Friedel oscillations 
and preserves the x~ 3 ^ 2 dependence of their amplitude 
far from the barrier (see Appendix B). That's why the lin- 
ear temperature dependence and change of sign of slope 
of G(T) should be robust with respect to the exact shape 
of the barrier. 

In summary, we calculated the conductance of a narrow 
and short quantum point contact at nonzero tempera- 
ture taking into account the electron-electron interaction. 
The conductance linearly depends on temperature and is 
proportional to the fourth power of the contact size, and 
the relative correction does not depend on the contact 
size. The sign of the linear temperature-dependent term 
depends on the competition between direct and exchange 
interaction. Measurements of the slope of its tempera- 
ture dependence allow one to determine the parameters 
of electron-electron interaction. 



unknown function into a Fourier integral 



dk 



2tt 



(Al) 



and then search a solution in the form of outgoing waves 
tpt(x,k y ) = c\ x e " x . The constant c\ is easily ex- 
pressed from the boundary condition 

00 

ci= / dy' e ik w' X {Q,y'). 



We sequentially substitute these formulas into (|A1[) and 
obtain 



tpt(x,y)= / dy' x(0, y') K(x, y, y') 



where the kernel is given by 



K(x,y,y') 



1 H[ l Xk^x* + {y-y>y) 



-— kx 
2 



yjx 2 + {y - y') 2 



(A2) 



(A3) 



The condition of continuity of the derivative of total 
wave function at the gap may be written as 



dj>o(x,y) 1 | dip t {x,y) 

lx=-0 



dx 



dx 



_ dip t (x,y) I 

x=-0 ~ Q x lx=+0' 

(A4) 



where 



ikx x 



We substitute the kernel (|A3[) into (IA4[) and obtain an 
integral equation 
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This integral equation appears in a problem of a two- 
dimensional diffraction by a narrow sliti 2 ^. We are in- 
teresting in the case of k » kp, so the limit kpa <C 1 
is equivalent to fa C 1. The solution of the integral 
equation in the limit of fa 1 is given by 2 ^ 



x(o,y') 



a 2 — y' 2 



(A6) 



Appendix A: Calculation of transmitted wave 
function 



We use the approach described by Sommerfeld— for 
the problem of two-dimensional diffraction and solve the 
system (J7J assuming a continuity of the derivative of the 
wave function ip = V'o + V't a t x = 0, where ipo is the wave 
function in the absence of the gap and ipt is the correction 
in the lowest order in a contact size. We expand the 



We substitute (lA6l) and (lA3l) into (lA2l) and obtain 



it \ n 7 1 2 1 1 f™> Hi(kr) 
ipt{x,y) = -k x ka mdj, • 

It is conveniently presented in the form 



, , , vk m 2 
Vt{x,y) = — 2Vk~ a 



dlv p i(q x x+q y y) 

2ir 



(A7) 



'"" , (A8) 



8 



In a similar way we solve the system (IT3l) for the 
Green's function of the Schrodinger equation g t and ob- 
tain 

8 nr rr 



Appendix B: Friedel oscillations of electron density 

The oscillating electron density is given by 

n(r,n) = ^ + (n)V-(r)) = ]T f(e a ) C(ri) ^(r), 

(Bl) 

where V'a are the eigenfunctions of the non-interacting 
Hamiltonian. In the lowest approximation in the contact 
size we neglect the distortion of the Friedel oscillations 
of electron density due to the gap and consider them as 
arising from a solid barrier. Because of the translational 
symmetry in the direction parallel to the barrier we may 
use instead of a momentum p and the eigenfunctions 



(B2) 



^ p (r) = e ip * y [e ip * x + r(p x )e- ip * x ] 



With the help of distribution function f(p), Eq. (|B1[) 
may be rewritten in the form 

<t,ti)= j j^mr v {*x)M*)- (B3) 

In the case of a sharp infinitely high barrier, the reflection 
coefficient r(p x ) = — 1 for any p x and the substitution of 
ip p (x,y) into (|F33[) gives us after simple rearrangements 

f d 2 p 

n(r,r!) = A J — cos[p y (y - y x )] 

x sm(p x x) sin^xi) (B4) 

Then we go to cylindrical coordinates (p, tp), expand the 
product of trigonometric functions into a sum of four 
terms and perform the integration over <p. The remaining 
integral over ip is a sum of four integrals of the form 



tt/2 



dip cos(acosip + bs'mip) = tt Jo(\/ a 2 + b 2 ). (B5) 



-it/2 



Calculations give the electron density in the form 



dp 



n(r,r 1 ) = / — f(p)p \J (p\/{x - x x ) 2 + (y - yi ) 2 ) 



-Mp^{x + Xl ) 2 + (y-yi) 2 ) 



Upon integrating by parts, one arrives at 

oo 



X p 



Ji{P\/(x - xi) 2 + (y - y x ) 2 ) 
V \x - xi) 2 + (y - yi) 2 

_ Ji(?V(x + .xi) 2 + {y-yi) 2 ) 
V(^ + a;i) 2 + (y-yi) 2 



(B6) 



In the case of ri = r the Friedel oscillations of density 
depend only on i, so 



n(r, r) = n{x) 



I 

2^ 



dp - 



df 
dp 

2 

P 

~2 



pJi(2px) 



2x 



(B7) 



In the case of a smooth yet impenetrable barrier the 
wave function ip p away from it still may be presented 
in the form (|B2[) with a momentum-dependent reflection 
coefficient r(p x ) = — exp[i5(p x )]i where 8{p x ) presents 
the phase shift of the reflected wave with respect to the 
case of zero boundary conditions at x = 0. We substi- 
tute (|B2I) into the Eq. (|B3|) and consider the coordinate- 
dependent oscillating part of electron density. With the 
help of cylindrical coordinates (p, ip) we can write it in 
the form 

dpf(p)p J dtp 



x cos(2px cos tp + 5 (p cos ip)). (B8) 

At large distances from the barrier 2px 3> 1 we can esti- 
mate n osc (r, r) using the stationary phase method 



2tt z 



dpf(p)p 



— cos(2px - 7r/4 + 5(p)) + O 
px 



V(2px) 3 / 2 



(B9) 



If dd(p)/dp <C 2x, it is possible to integrate this equation 
by parts and to obtain 



(r,r) 



1 



4(7ra) 3 /2 



dp 



dp 



x p 1/2 cos(2px + 5(p) +tt/4). (BIO) 



It is easily seen that at low temperatures and far from 
the contact, n osc exhibits the same asymptotic behavior 
as the oscillating part of Eq. (|B7[) except for the phase 
shift S(pf)- 



9 



2DEG 




FIG. 5: A model of experimental structure. 



Elucidate now the conditions for validity of (|B10[) . If 
the screening in the system is two-dimensional and the 
potential of the barrier falls off according to a power 
law, d5(p)/dp ~ Xo(p), where xo(p) is the classical 
turning point for the electrons at the barrier. Even if 
the width of the barrier 2xq(pf) is of the order of the 
Fermi wavelength, there is a large interval of distances 
Pp 1 < x < v F /T where Eq. (|B10I) holds, and it is 
precisely this interval that dominates the temperature- 
dependent contribution from Friedel oscillations to the 
conductance. 



to the in-plane coordinates and integrate it with respect 
to z over a small vicinity of z — to arrive at an equation 



dc/)(q, z), z=+0 , , 



47reo 



(C2) 



Here K2 — 2/ as is the inverse two-dimensional screening 
length. We assume that the layer of 2DEG is thin and 
the potential is continuous at z = 0, i.e. 

(f>(q,z = -0) = ^q,z = +Q). (C3) 

We solve Eqs. (|C2[) and (|C3j) with boundary conditions 

(C4) 
(C5) 

(C6) 



<f>(q,d) = Q 
0(9,-oo) = 
and write the potential in the form 



0(9,0) = 



1 



[coth(qd) + l]q + 2k 2 ' 



The potential in the coordinate space is obtained by in- 
verse Fourier transform of Eq. (|C6I) and may be conve- 
niently expressed in terms of a dimensionless coordinate 

x = qru 



Appendix C: Screened Coulomb potential 

We calculate Coulomb potential screened by a gate and 
two-dimensional electrons. Consider the system shown 
in Fig. [5] with a positively charged particle eo at point 
(0, 0, 0). The total potential induced by the charged par- 
ticle, the two-dimensional electrons and the gate satisfies 
the Poisson equation 



4-7T 

V 2 0(r || ,z) = -[p^ + p 1 ' 



(CI) 



eo<5(i"||)(5(z) is the density of the particle 
charge and p md — — e 2 V2S(z)(j)(r) is the density of the 
induced charge calculated in Thomas-Fermi approxima- 
tion. We take a Fourier transform of (|C1[) with respect 



r lh°) 



2e 



er\\ 



dx 



J (x) 



coih(xd/r\\) + 1 + 2k2T\\/x 



(C7) 



In experiment s 16 ! 17 the case of d 3> n^ 1 is realized. An 
evaluation of the integral (|C7[) gives for different limiting 

cases 



(r ,0) 



eo 



9 ^ ■ 
eg 

er||' 



r-| 1 > d » k 2 1 



d > ri > k, 



d > k 2 1 > ni 



(C8) 
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